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It is shown that an asymmetric spin-orbit interaction is appreciably renormalized by the effect 
of critical spin fluctuations. As an explicit example, the Rashba-type interaction is predicted to be 
suppressed by anti-ferromagnetic critical fluctuations near the hot line (or point), leading to a defor¬ 
mation of the Fermi surface near the hot line (or point) around which the band splitting is minimized, 
while ferromagnetic critical fluctuations are shown to enhance the Rashba-type interaction every¬ 
where on the Fermi surface, leading to an increase in the band splitting. It is also found that the 
many-body mass renormalization, which diverges towards the magnetic quantum critical point, is in 
the opposite sign for the two split bands owing to the Rashba-type interaction. These predictions are 
observable through the de Haas-van Alphen experiment in principle. 


The role of asymmetric spin-orbit (ASSO) interaction in giving rise to intriguing physical phe¬ 
nomena of strongly correlated electron systems has attracted much attention in the past decade. In 
particular, it has been discussed in relation to the superconductivity of non-centrosymmetric systems 
without the inversion centers of crystals in the field of heavy electron systems . 1,21 It is also an essential 
ingredient for the surface states of topological insulators, which have been extensively discussed for 
the past several years since their discovery by Kane and Mele. 3) 

It is an old problem how many-body effects, owing to the on-site Coulomb interaction, affect the 
atomic spin-orbit (SO) interaction in orbitally degenerate d-electron systems . 4,51 Recently, phenomena 
of the growth in the spin Hall effect have attracted much attention, and theoretical proposals have been 
published from the viewpoint of understanding that these phenomena are due to the effect of the growth 
of the atomic SO interaction . 6,71 

On the other hand, it has also been discussed how the ASSO interaction of Rashba type influences 
crucially the effect of critical spin fluctuations (CSFs) on the superconducting transition temperature 
near the antiferromagnetic quantum critical point . 8,91 Furthermore, it has recently been discussed how 
the Kondo effect is influenced if the conduction electrons are subject to strong ASSO interaction. 
For example, they are concerned with behaviors of magnetic impurities on the surface of topologi¬ 
cal insulators 10,111 and the Kondo effect in the system where conduction electrons are subject to the 
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Rashba-type spin-orbit interaction. 12,13) 

However, it has not been clarified how CSFs affect the ASSO interaction. This effect possibly gives 
rise to a deformation of the Fermi surface near the magnetic quantum critical point. In this Letter, we 
report on a perturbation calculation of a correction to the ASSO interaction under the influence of 
CSF modes, revealing the fact that CSFs, ferromagnetic (F) or antiferromagnetic (AF), give rise to 
a singular renormalization of the ASSO interaction. This fact suggests that the ASSO interaction is 
renormalized strongly by many-body effects. 

We discuss the renormalization of the ASSO interaction due to CSFs near the magnetic quantum 
critical point by a perturbation calculation. Such a renormalization is expected to be visibly large 
considering that CSFs considerably enhance nonmagnetic impurity scattering. 14) We demonstrate that 
the ASSO interaction is markedly influenced by CSFs, leading to the deformation of the Fermi surface 
(FS). We also note that the dynamical effect on the renormalization for ASSO interaction gives rise to 
a nontrivial mass renormalization of quasiparticles and the deformation of the FS. 

To provide an explicit discussion, we adopt a single-band model action, on the basis of the idea of 
an itinerant-localized duality model with tetragonal symmetry but without an inversion center of the 


crystal as follows: 15) 

A = Af + A s + A; n t + A so , (1) 

A f = - Z f^(-ie n )\G^(k, is n )\- ] f klT (is n ), (2) 

k,cr,n 

A s = - Z S- q (-;V m ) • S q (iv m )\xo(iv m )~ ] - 7(q)], (3) 

q,m 

^int = — ^ ^ ^ ^ j ./—k—q ,cr^~^ n ~ ^m) & cr,<r' f kcr' (i&n) * (4) 

k,q,<x,cr' n,m 

-3so = O' Z ZrT( k ) • fro-o-'fko-', (5) 

k,<x,cr' 


where f ktT (/J r ) is the annihilation (creation) operator for the quasiparticle, with the momentum k and 
the “spin” cr, renormalized by local correlation, and S (q) is the local spin density operator constructed 
from high-energy states of correlated electrons. Note that “spin” stands for pseudo-spin specifying a 
degeneracy of the Kramers doublet in general. Hereafter, we call it as simply spin. The Green function 
} (k, is n ) of the quasiparticle is given by 

G^(k, ie„) =* -— —r-, (6) 

sk 

where is the dispersion of the quasiparticle measured from the Fermi energy, s n = (2 n + I )nT is the 
fermionic Matsubara frequency, and z is the renormalization amplitude due to the effect of the local 
correlation, xo(jvm) (with v m = 2 rrmT being the bosonic Matsubara frequency) and 7(q) in Eq. (3) rep¬ 
resent the local dynamical spin susceptibility and the exchange interaction among localized parts of 
spin degrees of freedom, respectively. Aj nt [Eq. (4)] represents the exchange interaction between quasi- 
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particles and local spins. Finally, A so [Eq. (5)] represents the ASSO interaction: 1,21 a is the coupling 
constant of the ASSO interaction, and cf = (cr x , cry, cr-j is the vector spanned by the Pauli matrices. 
The vector y (k) = (yt^yt^yt) in Eq. (5) obeys y(-k) = -y(k), thereby breaking the inversion 
symmetry and causing the ASSO interaction. The main purpose of this study is to demonstrate how 
the ASSO coupling constant a is renormalized by CSFs. Hereafter, to provide explicit discussions, 
we assume that the ASSO interaction is the Rashba-type SO interaction given by 1 61 

y (k) = (sin k y a, - sin k x a, ()), (7) 


where a is the lattice constant in the tetragonal plane. However, it remains valid also for another type 
of ASSO interaction 17 - 181 that the ASSO coupling constant a is considerably influenced by CSFs. 

The lowest-order correction to the coupling constant a of the ASSO interaction due to CSFs is 
depicted by the Feynman diagrams shown in Fig. 1. The off-diagonal (in spin space) self-energy of 
quasiparticles due to this correction (k, ie n ) appearing in Fig. 1(a) is given by 


ElU (k - is n) = -^ T Yj a \ y ( k - k '> ten - ten') [G (0) (k\ lS n ')\ 2 , (8) 

k' n' 

where xu is a part of the spin-fluctuation propagator y s = (y^ -y^)/2. 

The propagator xu is expressed in terms ofy s and the charge-fluctuation propagator y c = 2(xjj + 
y-jq) asy-jq = ~Xs+(Xc/ 4), in general. Near the magnetic critical point, the propagatory-jq is dominated 
by the spin-fluctuation propagator y s and is given by 14,19-211 


au (q A<An) - -- 


(0) 

Xq 


n z -, (9) 

p + A (q - Q)“ + C q \a> m \ 

where a> m = 2m nT is the boson Matsubara frequency, rj is the distance from the magnetic critical 
point, Q is the magnetic ordering vector, andy^ 1 is the noninteracting static spin susceptibility at the 
magnetic ordering vector Q. Note that rj > 0 represents the paramagnetic state and // = 0 corresponds 
to the magnetic critical point. 20 - 211 

Substituting Eq. (9) into Eq. (8), we obtain 


FnTi ( k < te n ) = ^ (cr y ) n r y (k, is n ), 


( 10 ) 


with 

T v (k, is n ) = 




aAT LL^ k,) : 


A, 


(0) 


1 


k' n' 


ri + A(k - k' - Q)- + C| k - k '| |e„ - s n '\ (is n > - £ k ') 2 


(v = x,y), 


( 11 ) 


where A = zA. The other off-diagonal self-energy Enjq (k, is n ), shown in Fig. 1(b), is given by 

Injq (k, is n ) = ^ (oy)*^* (k, -ie „), (12) 

v=x,y 

where we have used a general relation 


G iT (k,/e„) = [G n (k, -is n )T, 


(13) 
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Fig. 1. Feynman diagrams for the lowest-order correction to the coupling constant a of the ASSO interaction 
due to CSFs. The wavy line represents the spin-fluctuation propagator;^, the solid line the Green function 
G of the quasiparticles, the filled circle the coupling constant A between spin fluctuations and quasiparticles, 
the upward arrow the up spin, the downward arrow the down spin, and the cross the ASSO interaction. 


as shown in Sect. 1 of Ref. 22. 

First, we discuss the effect of three-dimensional (3D) AF quantum critical fluctuations (QCFs) in 
which C q is ^-independent, i.e., C q = C. We also discuss the case where the external wave vector k is 
on the FS. To begin with, let us define k* as k* = k - Q and investigate the case where k* is on the 
spherical FS; i.e., k is on the hot line where (A,q = Note that the results obtained in the present 
study are valid for any shape of FS as long as Eq. (9) is valid, although the calculations become far 
more tedious, ffere we assume that the origin of the AF order is caused by a commensurate AF- 
exchange interaction through an incoherent process in the meaning of the itinerant-localized duality 
model so that the ordering vector Q is not influenced by QCFs , 15 - 23 ) 

Then, the frequency-independent terms in the denominator of Eq. (9) are transformed into 

ij + A(k - k' - Q) 2 = 11 + A(k* - k') 2 - 4r(£' 2 + 2 b? + a), (14) 

V F 

where stands for (A, and a and b are defined by 

Vp 22 

a = —T] + 2&pVp(l - cos 8), (15) 

jCX 
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and 


b = &fVf( 1 - cos 9), 


(16) 


where cos 9 = (k*-k')/(|k*||k'|), and ky and vp are the Fermi wavenumber and velocity of quasiparticles 
on the FS, respectively. Thus, T y (k, is n ) [Eq. (11)] is expressed as 


F y (k, is n ) = -aA 2 y v ( k-Q) T £ £ —^ J 


X, 


(0) 


1 


f ,:2 + 2 b?+a + C*\e n - s n >I (£' - ie n >Y 


(17) 


where C* = C/A. Since the singular behavior in T v appears through the contributions from the 
regions g - 0 and cos 9 - 1, in deriving Eq. (17), we have used the approximation k' w k* = k - Q; 
i.e., y v (k') ~ y v (k - Q). Similarly, the density of states of quasiparticles, is approximated by 

that at the Fermi level, Ny. Therefore, r v (k, is n ) is reduced to 


r v (k, is n ) = -Jy v { k - Q) 

where J is defined by 




i 


i 


2 + 2b + a + C*\s n - E n '\ {%' - is n ') 2 


J = a~A 2 N ¥X Q V -j- 


, (18) 


(19) 


Performing the integration with respect to in Eq. (18) with the use of the contour integration, 
we obtain 

d (cos 9) 


T v (k, is n ) = -7y v (k- Q) 


r, 


7TT 


y {- 


i 


2 

1 


1 


1 


n '>0 V 


( 20 ) 


Va+ - b 2 [b + i( y]ci+ - b 2 + e„')] 2 _ b 2 [b - i( yja- - b 2 + e„')] 2 

where a± = a + C*\ + s n - e„»| Here we note that Eq. (20) is valid regardless of the sign of e„, positive 
or negative, although we consider the case with s n > 0 hereafter. Equation (20) is rearranged and 
approximated as 

d( cos 6) 


T y (k, is n ) - T v (k, is n ) + Jy v (k - Q) 


nT 


! z 

I n'>0, n'±n 


i: 


i 


i 


i 


i 


+2ib L 

n'>0 


Vd7 (Vd++ £«') 2 yjdl + E n ') 2 

11 11 


( 21 ) 


yjdl (+ £n ') 3 s/d^ (^/^Z + s n ') 3 
where f v (k, is n ) is the contribution from n' = n > 0 in the first term in the brace of Eq. (20), which 
corresponds to the static term of fluctuation and should be treated separately, as discussed in Sect. 3 of 
Ref. 22. In deriving Eq. (21), we have neglected the terms of 0(b 2 ) because the singular contribution 
to T y arises from the regions cos 6 - I and b 2 oc (1 - cos 9) 2 . Then, the retarded function (k, e + id) 


5/12 






















J. Phys. Soc. Jpn. 


Letter 


at T —* 0 in the static limit (e —> 0) is given by lim^o f v (k, inT): 


r* (k. id) = (k, id) +/y v (k - Q) 


r 1 encase) r d£ , 

J-i 2 Jo 


Va + C*s'( Va + C*e' + £ ') 2 ’ 


( 22 ) 


where the term f ^ (k, id) is given by the first term in the bracket of Eq. (23) in Sect. 3 of Ref. 22. 
By changing the integration variable from e' to y = \la + C*e' / y[a, Eq. (22) is reduced to a more 
compact form as 


r* (k, id) = ar af (? 7 )y v (k - Q), 

where the excess contribution a a f to the SO coupling from the spin fluctuations is given by 


a Ad) = 


T* (k, id) 2J f l d(cos8) 1 
7 v(k - Q) 


+ 


u r d(cos 6) i r° 

oL 2 Jj 


dy 


1 






2 ' 


(23) 


(24) 


With the use of the expression for a [Eq. (15)], it is easier to perform first the integration with respect 
to cos 6 and then with respect to y. According to Sects. 2 and 3 of Ref. 22, the result for a a f up to 
0( yjrj) is given by 


a-Ad) 


2k 2 v 2 


Am* A \ 2 VA 


v F C 


C 


An + 


Ak l nT 
k F v F An 


+Q(n ) 


(25) 


where m* is the effective mass of quasiparticles and Fix) is defined by 


F(x) = log - 


1 


VT 


: log 


x+ i - VT 


+ X" 


+ X" 


X + 1 + Vl + X 2 


(26) 


The function F{x) is an increasing function and positive definite for x > 0, and F( I) ~ 0.55. Therefore, 
T^ (k, id) has a cusp singularity at n = 0 . 

The s n dependence of T v (k, is n ) for s n > 0, up to 0(e„), is given in Sect. 3 of Ref. 22 as 

T v (k, is n ) = Y* (k, id) + y v (k - Q) h (; is n ) + 0(4), (27) 


where h is defined by 


h = J 



nT 

4 + 


1 

4(k F v F ) 3 


log 



(28) 


We note that h is divergent because 77 ex 7' ,/2 just at the 3D-AFQCP. 24) In the case where the result 
of the first-order perturbation calculation is divergent, we need to consider the effect of higher-order 
corrections with respect to CSFs in general. This problem is beyond our scope and will be left as a 
future work. 

Substituting Eq. (27) into Eqs. (10) and (12) and performing the analytic continuation is n —* e + id, 
we obtain 


(k, e + id) - y (k - Q) • (<T) n [a af in) + hs ], 


(29) 
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and 

Ij iT (k, £ + iS) - y(k - Q) ■ (d-)^[ a ai (rj) + hs\. (30) 

Then, the retarded Green function matrix G R for quasiparticles, renormalized by QCFs, is given by 

G R (k, s) = [(£ -4)1- a 7 (k) ■ & - S R (k, e)] 1 , (31) 

with 


I R (k, e) = Sf (k, e) 1 + I R (k, e). (32) 

Here, S R (k, s) is the self-energy of the quasiparticles due to the effect of spin fluctuations without 
considering the Rashba-type SO interaction. 

According to Sect. 4 of Ref. 22, by diagonalizing the Green function matrix G R , the dispersion of 
quasiparticles is given by 

Ik = [4 + If (k, 0) ± |ay (k) + oraf (rj) 7 (k - Q)|] (33) 


where 



dlf (k, e) 
ds 


s=0 


h [ay( k) + a a f(q)y(k - Q)] ■ y(k - Q) 
\ay(k) + a a f(q)y(k - Q)| 


(34) 


It can be seen from Eqs. (33) and (34) that the band (0) with the energy = [4 + lf(k, 0)]2k, 
where 2k is equal to 2^ with h = 0, is split into the two bands (±) with the energies and the term 
of a s {q)y (k - Q) is newly involved into the coupling constant of the Rashba-type SO interaction due 
to CSFs. The FS is deformed by an effect of 24(k, 0) in general. Therefore, precisely speaking, the 
calculation of Z R (k, s) in Eq. (32) should have been performed with the use of quasiparticles with 
this deformed FS. However, we are interested in the effect of £f(k, e) on the band splitting due to 
the ASSO interaction, and we have taken an approximation that the spin-fluctuation propagator (9) is 
maintained intact. Note also that the divergences in h [Eq. (28)] can make either of fe) negative, 
which is unphysical. To remedy this fault, we need to consider the effect of higher-order corrections 
with respect to CSFs, as pointed out above. 

The derivations above are restricted to the case where k is on the hot line(s) on the FS without the 
ASSO interaction. To discuss a deformation of the FS due to the ASSO interaction renormalized by the 
AF-CSF, the effect of the distance of the position of k from the hot line on the FS is parameterized by 
shifting tj to rj + Aq} n in Eq. (25), 14) where q m is the perpendicular distance between k and the hot line 
on the FS without the ASSO interaction. Then, in the case where the ordering vector is commensurate, 
i.e., Q = (n/a,n/a,Tr/c), so that y v (k - Q) = -y v (k), the band energy ^ in Eq. (33) is approximated 
as 


£k - 


4 + if (k, 0) ± a - ff a f (q + Aq%) 11 y (k)| 


'k ’ 


(35) 


where 2^ is given by modifying h in Eq. (34) with the replacement of // by // + Aqy n in Eq. (28). There¬ 
fore, according to Eq. (34), we expect that the effect of dynamical mass renormalization is opposite in 
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the two split bands: the effective mass in band (+) is enhanced while that in band (-) is suppressed. 

The singular behavior given by Eq. (25) fades away at positions on the FS far from the hot lines. As 
a result, the suppression of the band splitting is confined to the region near the hot line, as shown in Fig. 
2(b). Thus, it means that the deformation of the FS is caused by the band splitting under the influence of 
critical AF spin fluctuations. In particular, the suppression of the band splitting seems to be followed 
by an AF-gap opening around the hot line(s), which is reasonable from the physical viewpoint. In 
the “AF ordered state”, the ordering vector Q itself should be determined self-consistently by taking 
into account the effect of the enhancement of the ASSO interaction a. Therefore, it is possible that a 
resultant ordering vector Q can be incommensurate even though the critical mode is commensurate in 
the paramagnetic state. This problem will be discussed elsewhere. 



Fig. 2. Fermi surfaces split by the Rashba-SO interaction for (a) the system without the effect of CSFs, (b) 
the system with critical AF spin fluctuations, and (c) the system with critical F spin fluctuations. The circle 
indicated by the dotted line represents the Fermi surface of the band (0), the large circle indicated by the 
solid line with the (±) represents the Fermi surface of the band (±), the long leftward arrow represents the 
AF-order vector Q, the short arrow represents the direction of the spin, and the small open circle indicated 
by the solid line represents the hot spot. 


There exists a singular contribution to the renormalization amplitude 2^ near the criticality (?/ = 0) 
through a singularity included in li [Eq. (28)]. Another nontrivial aspect suggested by Eqs. (33) and 
(34) is that the mass renormalization effects [through the factor ] for the two split bands are in 
the opposite sign, as explicitly shown in Eq. (34). This effect can be observed through the de Haas-van 
Alphen (dHvA) experiment in principle. 

Note that the aspects discussed above are expected to survive even if the ordering vector Q is not 
commensurate as long as the deviation from the commensurate one is not very large. 

Next, let us discuss the effect of 3D F-QCP in which C q in Eq. (9) is proportional to q~': C q = 
y/lk^C/q, with q - |k - k'| =* V2&f VT - cos 6 and C being constant. Then, after calculations similar 
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r, (k, id) = a ( (i])y v (k), 


where af is given by 


<ff(?7) = 


r*(k ,iS) 

7v(k) 


+2 J 


i: 


d(cos6) VT - cos 6 
2 C* V« 


J 1 


cly 


y + ^— i (y 2 


n 2 ’ 


1 ) 
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(36) 


(37) 


where C* = v 2 C/A. As shown in Sects. 5 and 6 of Ref. 22, by calculations similar to those deriving 
Eq. (25) from Eq. (24), the (k, id) is, up to a logarithmic accuracy in //, reduced to 


aftn) - J 


1 


4^pVp 


4 Vim* A 
C 




Ak 2 

F 


G 


4 V2 m*A 
C 


log 


16Ak 2 

F 


+ H 


2 Vim* A 
C 


log 


2 

e r/ 


l6Ak£ 

F 


Ak l JlT 
k ¥ v ¥ ^ 


where the functions G(x) and H(x) are defined by 

T 2 

G(x) = 


8(1 + x 2 ) 


l x 

1 + - + log - 
^ 2 


Fix) 


and 


H(x) = 


x^_ r 

~ J i 


cly ■ 


y 2 — i 


i)]- 


(38) 


(39) 


(40) 


[y + x(y 2 

Both the functions G(x) and H(x) are positive definite. 

The coefficient ( 77 ) increases sharply as the critical point is approached, i.e., // ^ 0. so that it 
exhibits a shaip cusp as a function of 77 . As a result, the band energy ^ in Eq. (33) is approximated as 


£k - [& + ( k > 0) ± \a + af(ii)\ |y (k)|] z£, 


(41) 


where is defined by the same expression as Eq. (34) with the same expression as h [Eq. (28)], as 
discussed in Sect. 6 of Ref. 22: 


h = J 


l Ak f nT 


+ 


1 


(^ F vp) 4 Vn 4(^ f vf) 3 


log 


4 Akl 
rje 


(42) 


This means that the coupling constant of the Rashba-SO interaction is enhanced due to the critical F 
spin fluctuations, resulting in an increase in the splitting of the band ( 0 ) into the bands (±), which we 
can easily see if we compare Fig. 2(c) with Fig. 2(a), in contrast to the AF critical fluctuations where 
the band splitting is suppressed. This marked difference arises from the fact that y v (k - Q) = y v (Q) 
for Q = 0 and y v (k - Q) = -y v (Q) for Q = ( n/a , n/a, n/c), because of Eq. (7). As in the case of 3D- 
AFQCP, we expect that the effect of dynamical mass renormalization is opposite in two split bands: 
the effective mass in band (+) is suppressed while that in band (-) is enhanced, which is in contrast 
with the case of 3D AF-QCR 

We note that the singular behaviors in Eqs. (25) and (38) become more prominent in two dimen- 
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sions (2D) than in the 3D case. These singularities are obtained by replacing the integration with 

r*2 71 

respect to cos 6 in Eqs. (24) and (37) by L dtp/In and (1 - cos 6) by (1 - cos ip) with <p being the 
angle between k* and k'. As shown in Sects. 7 and 9 of Ref. 22, in the case of the AF critical point we 
obtain 


a ai {ri + Aq 2 ) = 


-f 


2m* A 644 A 2 4 Akl 

log-!- + 


2k 2 Vp i nvpC q + Aq 2 k F vp q + Aq, 


(43) 


instead of Eq. (25), and as shown in Sects. 8 and 10 of Ref. 23; in the case of the FM critical point, 
instead of Eq. (38), we obtain 22 ^ 


a f (q) = 


2&pVp 


27 


y/2m*A\ 2V2A _ lAklr 


where the function I{x) is defined by 


I(x) = x 


J 1 


C 


dy 


nCv p 


V?7 + 


kpVp 77 


+ 0(? 7 ) 


1 


y + x(y 2 - 1 ) 


(44) 


(45) 


y 3 / 2 [y + 2 x(y 2 - l )] 3 / 2 

In the cases of both 2D-AFQCP and 2D-FQCR the coefficient /i, in the renormalization amplitude 2^ 
[Eq. (34)], is replaced by 


h = J 


Akl 


T 
— + 


Akl 


1 


(£fVf) 4 7 4(/:fVf) 3 y/q 


(46) 


as shown in Sects. 9 and 10 of Ref. 22. 

Namely, the extent of singularities in or a f, ay, and h is strengthened compared with the 3d case 
because q oc -T / log T in 2D-AFQCP 25 ^ and q oc -T log T in 2D-FQCP. 26) In particular, in the case 
of 2D-AFQCP, the correction to the Rashba-type coupling constant on the hot line (q m = 0) diverges 
logarithmically as the critical point is approached, which suggests that the band splitting due to the 
Rashba-type SO interaction becomes vanishingly small on the hot line. 

In conclusion, on the basis of perturbation in CSFs, we have shown that there exists the off- 
diagonal (in spin space) self-energy of quasiparticles, which gives rise to a considerable renormal¬ 
ization of the ASSO interaction. As a result, the ASSO has been shown to be decreased near the hot 
line(s) on the FS near the AF-QCP while it is uniformly increased over the whole FS near the F-QCP, 
leading to the deformation of the FS in different manners for AF-CSF and F-CSF. We have also found 
that the effective masses of quasiparticles in the two split bands are influenced by the ASSO interaction 
in opposite directions. These predictions can be observed through the dHvA experiment in principle. 
Higher-order effects in CSFs, and a self-consistent treatment of the ASSO interaction and CSFs, de¬ 
serve further investigations. This type of effect on the ASSO interaction arising from enhanced spin 
fluctuations may be possible in other physical situations, which are also interesting subjects to be 
explored in future studies. 
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1. Derivation of Eqs. (12) and (13) in the text 

Off-diagonal components of fermionic Matsubara Green function matrix in spin space, G(k, r), 
are defined by Gjq(k, r) = -(ak-f(T)a^(0)> and Gp|-(k,r) = -(akiCOfl^lO)), for 0 < r < {3. Then, it is 
easy to show that the following relation holds: 

G iT (k,T) = [G n (k,r)]*. (1) 

The Fourier transform of the off-diagonal components, Gjq(k, r) and G-jq(k, r), are defined by 

r* 

G iT (k,i£„)= I d7V £ " T G iT (k,T), (2) 

Jo 

and 

rP 

G n (k,ie„)= I drd ie ' ,T G n (k,r). (3) 

Jo 

With the use of Eq. ([Hi, it is easy to show that the following relation holds: 

G iT (k,i £ „) = [G n (k,-i £ „)]*. (4) 

This is nothing but the relation Eq. (13) in the text. 

Since the inverse matrix of G(k, ie„) _1 satisfies the same relation as Eq. (@]). the off-diagonal 
components of the matrix self-energy, X(k, if;,, ), satisfies the same relation as Eq. ©■ Then, considering 
the relation (cr v )^ = (cr v )L , Eq. (12) in the text is immediately derived. 

2. Derivation of Eq. (25) in the text 

Let us introduce S 3 dAF(??) representing double integration part in Eq. (24) in the text. Substituting 
the expression (15) in the text into the expression of S 3 dAi : ( ? /)- .SJdArJ/) is explicitly written as 

X I 1 1 /"»00 1 

4 , dy-- _ --J. (5) 

i ^ sJn + K(] -X) Ji ^ y jf l + K{l-x){y 2 - 1 ) + y] 

m 
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where x = cos 6, fj = v^tj/A and K = 2k 2 v^. Changing the integration variable from x to u = 
yjfj + K( I - x), Eq. © is reduced to 

. yjfi+lK 


S 3dAF(?7) - 


1 r oo 

— I Au I dy- 

J l ’ U(v 2 - 


[£(y 2 -l)+y]' 


Then, the integration with respect to u is easily performed: 


S 3dAF(?7) - 


c* r° 

~~~K Ji 


dy 


1 


y 2 -l 


1 


where Bn and B\ are defined as follows: 


Bn = 


B i = 

Then, the integrand in Eq. ([7]) is transformed into 

Bo 


B\iy 2 - 1) + y Bo(y 2 - 1) + y 

V? 

C* 

yfn + 2K 


C* 


s 3dAF(V> - 


c* r a 

~~k J i 


dy 


B i 


There holds the following definite integral: 


r°° , b i 

J i dy yB(y 2 - 


y Bo(y 2 - 1) + y yfii(y 2 -l) + y 

1 2B + 1 - Vl + 4B 2 ' 

= - log - , log - ; 

Vl +45 2 2B+1+ Vl + 4B 2 , 


1 


i) + y 2 

With the use of definition (26) in the text, this definite integral is expressed as 

f dv ---= ~F(2B). 

J l y B(y 2 - 1) + y 2 


Therefore, Eq. (ITOl) is expressed as 


C* 


s 3dA F (?7) = [F(2flo) - F(2B[)]. 

In the limit of rj —> 0, Bo —> 0 and B\ —* VtF/C* = B*. Then, 

lim S 3d AF(77) = ~ [F(0) - F(2B*)]. 

t ]—>0 ZA 

It is easily shown that B(0) = 0, and B* = 2m*A/C since K = 2k 2 v 2 and C* = v 2 C/A. Therefore, 

C* I Am*A 

lim S 3d AF(?7) = —FT- 

2K \ C 

This gives the first term in the bracket of Eq. (25) in the text. The second term arises from F(2Bo) but 
F{2B\) gives only the term of 0{rj) except for the first term of 0(rj {] ). Namely, 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


B(2B 0 ) = log B 0 - 


1 


V 


l+«*„ 2 


: log 


2#o + 1 — - 

\jl + 4®i? 

2Bq + 1 + - 

v/l+ 4B o 2 


- - log 


l-Bo 
1+B 0 


- 2B n . 


(16) 


2120 ] 
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Then, substituting the definition © into Eq. (IT3T) . we obtain 

£_ 

2 K 

_ VF 

K K^A 

and the second term in the bracket of Eq. (25) in the text. 


5 3dAFfa) - 5 3 dAF(0) = -^E(2S 0 ) - -^-2fi 0 


y/rj + 0(rj), 


(17) 


3. Derivation of Eqs. (27) and (28) in the text 

Here we derive Eqs. (27) and (28) in the text for 3D-AFQCP. Neglecting the b 2 terms in Eq. (20) 
in the text, T v (k, ie n ) is expressed as 

d(cos 8) 


r v (k, ie n ) = Jy v (k - Q) 


i: 


’ T Z 

n'> 0 


l 


l 


+ 


l 


l 


+2i b 


Vd+ (V^+ + £«') 2 yfdZ {yJdZ + E n ') 2 

11 11 


yjal (y/o7 + e M ') 3 a/o~ ( y[aZ + s„>) 3 


( 18 ) 


where a ± = a + C*| ± s n - s n > \ with a given by Eq. (15) in the text as a ~ ( v 2 rj/A ) + 2k 2 v 2 (l - cos 8), 
and b = kyvy( I - cos 8), Eq. (16) in the text. Here, we have neglected terms of 0(b 2 ) because the 
singular contribution to T v (k, ie n ) arises from the region cos 8 ~ 1. We also note that the expression 
m is valid irrespective of the sign of s n , positive or negative. 

Since we discuss the retarded function of low-energy quasiparticles, we need a linear dependence 
of r v (k, ie n ) for s n > 0. In the limit T —* 0, the summation over s n > is approximated by integration 
with respect to s' as in Eq. (21) in the text, except for the term e„> = e n which corresponds to taking 
static limit ( a> m = 0) of spin-fluctuation propagator given by Eq. (9) in the text. The latter term is 
usually neglected since it usually gives no contribution in the limit T —> 0. However, it turns out that 
the most singular contribution to the s n dependence arises as shown below. 

Let us define the term arising from e n > = s n in the first and the third terms in the brace of Eq. (IT8l) 
(or Eq. (18) in the text) as f v (k, ie n ). Then, 


F v (k, ie n ) = /y v (k - Q) tiT 


i: 


d(cos 8) 


1 


1 


(yfa + £„)" 


+ lib 


1 


1 


yfa (yfa + e „) 3 


The expression in the bracket is expanded in e n up to linear order as follows: 

12 b 3 

“T77 - -;S n + 2l— - 2l b-rjrSn + 0(s n ) 


(19) 


( 20 ) 


With the use of a = i) + Zf(l - jc) and b = y/K/2(l - x), with x = cos 8, fj = (v^rj/A), and K = 2k 2 v£, 
Eq. (fl9l ) is reduced to 

1 2 E n 


T v (k, ie n ) = Jy v ( k - Q) nT 


f 1 dx 

J J 


+2i 


[fj + K(l -x)] 3 / 2 [rj + K{l -x)] 2 
(1 -x) „ [K (1 - x)s n 


2 [fj + K(l -x)Y 


6i 


2 [rj + K(l -x)] 5 ' 2 


( 21 ) 


?>m 
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The integrations with respect to * in Eq. (l2Tb are easily performed to give 


T v (k, ie n ) - Jj v ( k -Q)^(-^=- %+ 1 

* \ V^7 V 


j 2K 4i s n 

y/2K 8 *7 V2£ V>7 


( 22 ) 


where we have retained the most singular term in each contribution in the curly bracket. Substituting 
77 = Vp 77 /A and K = 2k 2 v^ into Eq. l!22T) . we obtain 

I A 1-2 . ^ 

i nT 


( 


r v (k,ie n )-/y v (k-Q) 


Ak l nT 


Akl nT 


2(kpVp ) 3 y/^ 2(kpVp ) 4 ?7 4(kpvp ) 3 


4Ak 2 

log-£ - ic, 


'< nT 


" (k F vp) 4 yfrj 


(23) 

Since 77 oc T 3 ' /2 just at the 3d-AFQCP , 11 the first and the fourth terms in Eq. (1U3T) vanishes as T ]/4 while 
the second term diverges as 7’~ l/2 . Although the first term in Eq. (1231) vanishes as T l; l in the limit 
T —7 0 at the criticality, this T dependence gives a sharper cusp than the second term of Eq. (25) in 
the text because -\J :77 oc T i/4 there. 

The singular behavior as Eq. (I23T) . which arises from the term s n > = s n in the first and the third 
terms in the brace of Eq. (ITSl) . does not appear for the self-energy I](k, ie„) in Eq. (32) in the text. 
This is because the denominators [Z 7 ± i( -\Ja± - b 2 + s n >)\ 2 in Eq. (20) in the text should be replaced 
by [b ± i( yja ± -b 2 + e„/)] for I[(k, i s n ) so that ( \[a + s n ) 2 and ( \[u + e n ) 3 are replaced by ( \fa + e n ) 1 
and ( \fa + £„) 2 , respectively, in the denominator of Eq. (fT9l ). so that the contribution vanishes in the 
limit T -> 0. 

The summation in Eq. (fTST) over s n > (except for the term s n ’ = s n in the first and the third terms 
in the brace ) can be approximated, in the region T —» 0, by integration with respect to s’ as follows: 


T v (k, ie n ) - T v (k, ie n ) = Jy v ( k - Q) 


i: 


d(cos 9) 


where L(s n , cos 9) and M(s n , cos 9) are given by 

~ nT d s' 


L(s n , cos 9) = 


-r 

Jo 


[L(s n , cos 6) + 2 iZ 7 M(s n , cos 0)], (24) 


1 


+ 


+ 


2 -\Ja + C*(s n - e')( + C*(s„ - s') + s') 2 

d Is' 1 

:„+ 7 iT 2 sja + C*(s' - e„)( sja + C*(s' - s„) + s') 2 

1 


f 

OS„ 

r°° As/ 

Jo — 


yja + C*(s' + s n )( sja + C*(s' + s n ) + e ') 2 


(25) 


and 


M(s n , cos 6) - 


rs„-nT 

JO 

Je n 

-f 


As' 


1 


2 yja + C*(s„ - s')( -\Ja + C*(s n - s') + s') 3 
50 d£_ 1 _ 

+nT 2 + C* (s' - £„)( - \Ja + C*(s' - s n ) + s’) 3 

As' 1 


2 + C*(s' + £„)( V« + C*{s' + £„) + s') 3 


(26) 


In deriving Eqs. (1251) and (l26l) from Eq. (fTST) . we have used the relation between a definite integral 
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with respect to s' and a summation over n' in the region T —> 0 such as 
nT ^ F{s n r,s n ,T) + nT^F{s n >\s n ,T) 


0 <n'<n 


n'>n 


\ ns n -nT 

2 J d s'F(e'-s n ,T) + 


- r 

2 Js, 


d e'F(e';E n ,T) 


s n +nT 


C7* 

1 + O | — T-pr 2nT e, 

a ilL 


and 


i r° 

7rJ V F(e„/;e„,r) w - de'F(e';e„, T) 

,rr 0 2 Jo 


1+<9 I 


(27) 


(28) 


where F{s n r,s n , T ) stands for the functional forms in the brace of Eq. (fl8l> . The explicit T depen¬ 
dence represents that of rj included in the parameter a = (Vp tj/A) + 2k 2 v 2 ( I - cos 8). In deriving the 
approximate relation in Eqs. (1271) and (1281) . we have used the trapezoidal rule and the estimation of its 
relative accuracy. 

By changing the integration variable from s' to s" = s n - s' in the first term in Eq. (l25l) . from s' 
to s" = s' - s n in the second term in Eq. (l25l) . and from s' to s" = e„ + e' in the third term in Eq. i 
L{s n , cos 8) is reduced to 

e„ d£ // ] 


L{s n , cos 8) = 


X 


nT 2 + C*E"( Vfl + C*S" + £,, - E") 2 

00 d e " 1 


r 

JnT 
r*°° A 

L 


2 Vn + C*s"( Va + C*e" + e" + E n ) 2 
d e" 1 


JnT 2 Va + C*e"( Va + C*s" + e" - s„) 2 

r £n ds^ _ i _ 

JnT 2 Va + C*e"( yja + C*e" + e" - s n ) 2 


(29) 


It is easy to see that the sum of the second and the third terms of Eq. (1291 ) are even functions in e„. 
The summation of the first and the fourth terms in Eq. (l29l ) is calculated in the following way. 

Let us introduce A X(s n ) representing the summation of the first and the fourth terms of Eq. (1291 ) . 
Then, it is reduced to 


A X(e h ) = 


f e " d s" 1 

1 

1 

JnT 2 Va + C*e" 

.(Va + C*e" + s n - e") 2 

(Va + C*s" — £,, + s") 2 


= -4 


X 


,e " d£' 

nT 


Sn ~ S 


2 [a + C*s" - (s n - £") 2 ] 2 ' 


By changing the integration variable from e" to v = e" / s n , we obtain 

l -y 


A X(s n ) = -2 e\ 


if 

JnT/ 


d y- 


(30) 

(31) 

(32) 


JnT/s„ ' [a + C*s n y - e 2 (l - v) 2 ] 2 
Here, we note that nT/s n , the lower limit of integration in Eq. (f32l) . must satisfy nT/s n < 1 because 
we are discussing the analytic continuation from the upper half plane of i s n . Then, up to the order of 
(9(e 2 ), the denominator of Eq. (1321) is approximated by a 2 . Therefore, A X(s n ) is given by 

„2 r 


A X(s„) = - 

a z 


i 2 nT {nT) 2 


5{20] 
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s 2 2 (nT ) 2 

=-T + nT —s n - 

a z a z 


a 


E n 2 

~js n . 

a- a z 


(33) 


where the term (nT) 2 /cr has been discarded because it gives only a negligible contribution to Eq. (fT8l) 
in the asymptotic region T ~ 0 as far as terms up to 0(s n ) are concerned. Therefore, up to 0(s n ), 
L(s n , cos 0) is given by 


L(e n , cos 6) 


£ 


ds‘ 


1 


nT Va + C*s"( Va + C*e" + s") 2 a 


+ nT —E n + 0(s 2 ), 


(34) 


where the first term is the contribution from the second and the third terms of Eq. (l29l) . 
Similarly, the expression Eq. (1261) is reduced to 

f £ » dfT _1_ 

JnT 2 yja + C*s"( Va + C*s" + E n ~ e ") 3 


M(e„, cos 9) = 


L 


°° d e " _1_ 

2 Va + C*e"( Vn + C*s" + s" + b „) 3 


JnT 

r°° d £ " 
Jn', 


1 


£ 


nT 2 y] a + C*S"( + C*S" + s" - E n ) 3 

" d e " 1 


I o ,- ,----• (35) 

JnT 2 yja + C*s"( Va + C*s" + E" - E n ) 3 
The summation of the first and the fourth terms in Eq. (1351) is calculated in parallel to the derivation of 
Eq. (1331) . Let us introduce A Y(e„) representing the summation of the first and the fourth terms of Eq. 
(I35T) . Then, it is reduced to 


A Y(s n ) = 


■n)= r 

£ 


1 


■ £ " d e" 

nT 2 Va + C*E‘ 


1 


1 


+ 


(Va + C*e" + En - e"? ( Va + C*e" - En + e") 3 

d s" 2 (a + C*s") + 6(g„ - e") 2 


J nT 2 [a + C*e" - (E n - s") 2 Y 
By changing the integration variable from e" to y = s" /e„, we obtain 

a + C*s n y + 3e“(l -y) 2 


A Y(s„) = 


n ) — J 

JnT / 


dy- 


(36) 


(37) 


hT/s n ' [a + C*s n y - s 2 ( 1 - y) 2 ] 3 

Then, up to the order of 0(s n ), the numerator and the denominator of Eq. (1771) are approximated by a 
and a 3 , respectively. Therefore, A Y(s n ) is given by 


A Y(s n ) = ^§(l - — 
a 1 \ s n 

E n nT 

a 2 a 2 

Therefore, by expanding the second and the third terms of Eq. 
given by 


(38) 

up to 0(s n ), M(s n , cos 0) is 


M(e„, cos 9) 


+ Sn _ r°° _ 3fn _ 

a 2 a 2 JnT Va + C*s"( Va + C*s" + s'') 4 


+ 0(s 2 ). 


6120 ] 


(39) 
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By changing the integration variable from s" to p = s"/a, the third term on the r.h.s. of Eq. (1391) is 
reduced to 


i r 
a^ir 


dp- 


3e„ 


(40) 


hT/a y/l + C*p ( y/l + C*p + yfcipY 
Singularity arising from this term is much smaller than that from the second term on the r.h.s. of Eq. 
(Bdl ) because the singularity stems from the vanishing behavior in a = (v^ij/A) + 2k 2 : v 2 : ( I - cos 8) near 
cos 8 - I and rj - 0. Therefore, neglecting the third term on the r.h.s. of Eq. (l39l) . we obtain 


M(e„, cos 8) =* + % + 0(s 2 n ). 


(41) 


Collecting the results of Eqs. (1241) . (l34l) . and (14TT) . we obtain the f v , (k, ie n ) up to 0(s n ) as follows: 


T v (k, ie n ) - T v (k, ie n ) + Jy v ( k - Q) f d(< ^ S ^ f ds" 

J- 1 1 YJitT 


1 


Va + C*£"( Va + C*s" + s") 2 

+ 0(s 2 n ), (42) 


„ 2 l nT e n 

+nT —s n + 2i b\ -— H—r 

a z \ a~ 

where a = (v 2 r]/A) + lk 2 v 2 {\ - cos 8) and b = k\:Vy.( I - cos 8). Then, the linear terms of T v (k, ie n ) in 
e„ is obtained as follows: 


T v (k, ie n ) - T v (k, ie n ) 


e„=0 


r v (k,i£ n ) - T v (k, ie n ) 


e„=0 


+Jy v (k - Q) 


x: 


d(cos 8 ) 


2 b 
nT —s n + 2—(l £„) 


■ (43) 


We note that, in the limit T —> 0, the first term in the square bracket of Eq. (|421 ) is reduced to the 
second term of Eq. (22) in the text. Integrations with respect to cos 8 on the r.h.s. of Eq. (1421) are 
easily performed: 


x: 


1 d(cos 0 ) 1 


Akl 


a 2 4(kpVp ) 4 


1 


1 


Akl 1 


t) j) + 4Akl J 4(kpvp ) 4 rj’ 


(44) 


and 


x: 


d(cos 8) b 


1 


log 


4Akl + rj 4Akl 


1 


4 Akl 

log —(45) 


-i 2 a 2 8 (kpvp ) 3 ( e r\ 4 Ak 2 + 77 J 8(/:fVf) 3 ' rje 

where we have retained the most singular terms and within the logarithmic accuracy to obtain the 
last approximate expressions. 

Finally, substituting Eqs. (I23T) . (1441 ) . and (1451 ) into Eq. (I43T) . we obtain 


r y (k, i£„) - T v (k, i£ n ) | £n= o + Jy v ( k - Q) 


'< nT 


+ 


1 


(kpvp ) 4 yfrj 4(kpvp ) 3 


log 


4 Akl 

rje 


(i E n ) + 0{sl). (46) 


Here we note that the term obtained from the first term in the square bracket of Eq. (1431 ) cancels out 
the second term in the bracket of Eq. (I23T ). so that the term proportional to s n disappears but only that 
proportional to \s n remains. This result verifies the fact that the first line in the brace of Eq. (fl8l) . i.e., 
its real part, is an even function in e„. Equation (1461 gives Eqs. (27) and (28) in the text. We also note 
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4. Derivation of Eqs. (33) and (34) in the text 


With the use of Eqs. (10) and (12) in the text, an explicit form of Eq. (31) in the text is given by 


[(5 R (k, e + id)] ‘ 


(47) 


e - & - Sf (k, s) -ay (k) • <x n - X R n (k, s) 

-ay (k) • <x iT - I R T (k, s) s-£ k - (k, e) 

where X R and are defined by Eqs. (29) and (30) in the text, respectively. For concise presentation, 
let us introduce the following notation: 

f v (k, e) = ay v (k) + y v (k - Q) [a af (if) + (, h + ig)e]. (48) 

Although g = 0 as discussed in Sect. 3 (e.g., as derived from Eq. (l46l) ). we here give a general 
expression for the dispersion of quasiparticles assuming g ± 0 in general. Then, Eq. (1471) is written in 
a concise form as follows: 


[G R (k, s + id)] ‘ 


e - & - 2 R (k, s) - Yj V =x,y fy(k, e)cr. 


Ti ) 


Zv=^ r;(k, eXcrJ*)* e — & — Sf (k, s) 


Dispersion of quasiparticles is determined by the condition that |G R (k, s + id)| 1 = 0: 


[e ~ & ~ (k, s)f = J] f v (k, s)crl l x £ [f(k, e)]>: 


4T 


(49) 


(50) 


v=x,y 


v=x,y 


The r.h.s. of Eq. (1501) is transformed into 

2 fv<rl l x £ f= (if/ + |f/) + i (f x f *) z , 

y=x,y y=x,y 

where we have abbreviated f v (k, s) to f v for conciseness. With the use of Eq. 
two terms on the r.h.s. of Eq. (I5TI) : 

,2 , if. |2 _ r_ /i,\ , ../I, /-»m2 


(51) 

we calculate the 


ir.vf + |r v |- = [ary (k) + ar a fy (k - Q)f + 2/iey (k - Q) • [ay (k) + or af y (k - Q)] + 0(s z ), (52) 


and 

i (f x f *) = lags [y (k) x y (k - Q)] z . 

Then, the condition Eq. (I50l) . determining the dispersion of quasiparticles, is reduced to 

<ffif(k,e) 


(53) 


e-&-Xf(k,0) 


ds 


£=0 


= ± {[ay (k) + ar a fy (k - Q)] 2 + 2 age [y (k) x y (k - Q)] z 
+2 hs y (k - Q) • [ay (k) + a at y (k - Q)]} 1 /2 . 
Therefore, the dispersion of the quasiparticles is given by 


(54) 


el - 


dlf (k, e) 


ds 


£=0 


h [ny(k) + or af y(k - Q)] • y(k - Q) | 
|ny(k) + or af y(k - Q)| 


8 ( 20 ] 
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= £k + (k, 0) ± |ay (k) + a af y (k - Q)| ± 


a{ge) [y (k) x y (k - Q)] 


I ay (k) + aafy (k - Q)| 
This gives the dispersion relation presented as Eqs. (33) and (34) in the text if g = 0. 


(55) 


5. Derivation of Eq. (38) in the text 


Let us introduce S idFiv) representing double integration part in Eq. (38) in the text. Substituting 
the expression (15) in the text into the expression of S 3 dF(? 7 )> S 3 di (^7) is explicitly written as 

Vl -x 


S 3dp(jl) = 


-jc 


dx 


2 c* y/fj + K(l - X) 


r o 

J i 


dy 


1 


-^y/rj + K (! - x)(y 2 - l) + y 


2 ’ 


(56) 


where x = cos (9, /'/ = v^r]/A and K = 2k 2 v 2 . Changing the integration variable from x to s = V1 - x, 
Eq. (1561 ) is reduced to 


*5 3dF(?7) = 


-I 


VI 


ds 


C* yjfj + Ks 


r 


dv- 


^PV-1)+? 


i2' 


First, we calculate S 3dF(0) which is given by 


S 3 dF(0) = 


= f 

Jo c* Vk J i 


dy 


Vk. 


—(y 2 - 1) +y 


2 ‘ 


Changing the integration variable from s to 7 = s 2 , Eq. (1581) is reduced to 


S 3dF(0) = 


2 C 


v r* r 

" VK Jo J 1 


dy 


1 


# z(y 2 -l)+y 


2 • 


Then the integration with respect to z is easily performed: 

1 


S 3dF(0) = 


_L r 
~2KJ, 

1 r 
~2KJi 


dy 


dy 


y2 ~ l [^f(y 2 -i) + y y 

2 y[K 
C* 


^(P-U+y 


Then, with the use of Eq. (fT2l) . we obtain 


S 3d F (0) = —F 


1 (4VF|1 1 (4x[2m*A 


4 K 1 C 


4 K F 


C 


(57) 


(58) 


(59) 


(60) 


( 61 ) 


where we have used definitions, C* = Cv^/A and K = 2k 2 v 2 . 

Next, let us discuss the lowest-order correction in 77 , which will turn out to be of the order of 
0(77 log 77 ). To this end, we calculate S} L \\ : (r]) - .S’3di40): With the use of the expression Eq. dSTb . we 
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obtain 


S 3dF(?7) - S 3dF(0) = “77 

C* 


I 


V2 


ds ■ 


yjfj + Ks'■ 


fSHi; 


r 


dv- 


syJfj+Ks 7 


(y 2 -l)+y 


d y 


1 


Mil A, 2 


(y -1) + y 


r\s± 

Jo c* 


yjfj + Ks 2 VK 


r 


+ r 

Jo c* Vk J i 


dy 


dy 


1 


i^ij+Ks 7 

C* 


(y 2 - i) +y 


l 


sy/rj+Ks 2 


i2 


Vk. 


12u,2 


(y 1 -l)+y 


.(62) 


c , -(y - i) +y 

It is easy to see that the second term of Eq. (f62l) can be expanded into the Taylor series so that this 
term gives only 0(t]) at most and can be safely neglected up to the logarithmic accuracy. Let us define 
the integration with respect y in the first term of Eq. (I62T) as L(s): i.e., 

1 


r 


L(s) = dy - 


( y/fj+Ks 7 


(y 2 - l)+y 


(63) 


Then, the first term of Eq. (1621) . denoted by X{rj), is expressed as 


X(TJ) = 


=±r 

C* Jo 


V2 


ds 


L(s). 


(64) 


yjfj + Ks 2 VK 

Changing the integration variable from s to z = s 2 , Eq. (l64l ) is transformed through partial integration 
with respect to z as follows: 


1 p 

XijD = - dz 

2C* VK Jo 

1 


\z + ± 


L{yfz) 


2 C*VK 


2 i 2+ !)- 2 + i! ioe 


fi+i- V2 


f 


dz 


z|z+ fH 4 f log 



yjz+ | + yfz ^ 

In the limit of rj —» 0, the first term in the brace of Eq. (I65T) . denoted by 2 yf~KX { 1 Ur]), is given by 


(65) 


x<,, W‘^ log H x , 1 To I < v5) - 


(66) 
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where, according to Eq. (1631) . lim^_>o L( V2) is given by 


lim L( V 2 ) = 
>7 - >0 J 1 


dy 




n 2 ' 


(67) 


After elementary integrations, an explicit form of Eq. (1671) is 


r 2 B* 

lim L( v2) =-— 

?/-»o l+4 £* 2 


1 + — +logB* -F(2B*) 


where 


2 B 


B* = 


= —G(2B*), 
B* 


2 VK 


c* 


( 68 ) 


(69) 


The last equality in Eq. (1681) is obtained using the function G(x ) defined by Eq. (40) in the text. 

The second term in the brace of Eq. (l65l) . denoted by 2 xfKX^Gj), is estimated as follows: Chang¬ 
ing the integration variable from z to w = Kz/fj, X <2] (ri) is expressed as 


A ( 2 ) (77) = - 


2 K 3 ' 2 


r2K/fj 

I dw 

Jo 


y/w(w + 1 ) - w + — log 


Vw + 1 - y[w 


Vw + 1 + 


dL( yjfjw/K) 


dw 


(70) 


Since the upper limit of the integration with respect to w in Eq. (I70l) diverges in the limit 77 —> 0, 
there is a possibility that a singularity in 77 arises. Therefore, the behavior of dL( \Jfjw/ K )/dw at w = 
2 K/f} —* 00 should be estimated. With the use of expression for L(s), Eq. (I63T) . we obtain 
dL( yJfjw/K) j 


dw 


c*Vk 


r 


dy 


y 


1 


fjy/w(w+l), 2 _ 1) + y 

C* VK J 


n3 


w + 1 


+ 


w 


w 


w + 1 


(71) 


Near the upper limit of w = 2 K/f} in the limit of // —> 0 , dL( \Jfjw/K)/dw is estimated as 


dL( y/ifw/K) 


2fj 


dw 


c*Vk 


r 


dy 


^(y 2 -D + y 


(72) 


With the use of definition (41) in the text, Eq. (1721) is expressed as 

dL( yjfjw/K) 2fj 2 


In the limit of w 


and 


dw c* x[K B* 2 

j, the following asymptotic formulas hold: 

Vw + 1 - yjw 


H{B ). 


log 


Vw + 1 + yfw 


- ~ log |4w|, 


yjwiw + 1 ) “ w - 


(73) 


(74) 


(75) 


Therefore, X (2 \rf), Eq. (170l) . is given by 
A ( 2 ) (77) 


2ri 


2K 3 / 2 c* VIB* 2 


H{B*) 


^2 K/fj 


dw(-i - ~ log |4w| 


Lj> ]og ‘-A x ±H(B-). 

4 K 3 / 2 S 8 K B* 


(76) 


nm 
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In deriving the last approximate result in Eq. (l76l) . we have discarded the term of the order of 
<9(p 2 log 77 ), and used the definition of B*, Eq. (l69l ) . 

Collecting the above results, Eqs. (1661) . (I 68 T) . and (1761) . near rj ~ 0, X(rj), Eq. (f65l) . is given by 

4 


1 77 


G(2fi*)logg+ //(£*)log ^ 


(77) 


4 /C3/2 c*B* 

With the use of definitions of B*, Eq. (l69l) . C* = Cvp/A, and fj = v^tj/A, Eq. (1771) is transformed into 
a compact form as 


X(rj) — 


_L^_ 

4 KAk* 


G 


4 V 2m* A 
C 


e 77 

log- — + H 


16 Ak^ 


2 y/2m*A 
C 


log 


2 

e rj 
16 Ak^ 


Then, S 3 dF(? 7 ) - S 3 dF( 0 ), Eq. (1621) . is given by 


S 3dF(?7) - S 3dF(0) 


_Ljl 

4 KAIc* 


G 


4 'Jim* A 


C 


eri 

log- + H 

16A& 2 


2 Jim* A 


C 


log 


2 

e ij 


\6Ak p 


(78) 


(79) 


Finally, adding Eq. (RTFl) and Eq. ( I79l ). the second term of Eq. (37) in the text is reduced to the first 
and second terms in the brace of Eq. (38) in the text. 


6. Derivation of Eq. (42) in the text 


Here we derive Eq. (42) in the text for the case of 3D-FQCP. We start with Eq. (fJ~ 8 l) . the same 
equation as that in the case of 3D-AFQCP except that C*, in the parameter a± = a + C*\ ± s n - s n >\, is 
not constant but C* = €*/ V1 - cos 8 with C* being a constant. Then, we obtain the same expression 
for f v (k, i£ n ) as Eq. (IT9l) . However, since Eq. (fl9l) does not include the coefficient C*, the expression 
Eq. m f° r f v (k, ie n ) can be also applied to the present case of 3D-FQCP: 


Ty (k, ie n ) = 7y v (k) nT 


r 1 dx 

J-i 2 


2 S n 


[fj + K(l - x)f! 2 [fj + K{l-x)] 2 


+ 2 i 


(1 - x) 


6 i 


2 [rj + K(l -x)] 2 V 2 [fj + K(l - x )] 5 / 2 
Then, the expression Eq. (1231) for f v (k, ie n ) is also valid in the case of 3d-FQCP: 


T v (k, ie n ) - 7y v (k) 


Ak l nT 


Akl 


nT 


inT 


2(&fVf) 3 yfn 2(^fVf) 4 7 4(/:fVf) 3 


log 


4 Akl 


1 Sn 


'K nT 


(k F v F ) 4 xfrj 


(80) 


(81) 


Since 77 oc T lA ' just at the 3D-FQCP, r) the first and fourth terms in Eq. (| 8 IT) vanishes as T 1 /3 while 
the second term diverges as T -1 / 3 _ Although the first term in Eq. (I 8 TT) vanishes as T l; ' in the limit 
T —> 0 at the criticality, this T dependence gives a sharper cusp than the second term of Eq. (39) in 
the text because 77 log 77 oc T 1/3 log T there. 

The summation over s n > (except for the term s n > = s n ) in Eq. (fl 8 l) can be approximated, in the 
region T —> 0, by integration with respect to s’ as discussed in Sect. 3. Since the final results for this 
contribution in the case of 3D-AFQCP, i.e., Eqs. (1431) ~ (l45l) . also do not include the parameter C*, 
they are also applied to the case of 3D-FQCP. 
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Then, T v (k, ie n ) is given by the same expression as Eq. (l46l) : 


T v (k, ie n ) - T v (k, ie n ) | e „ =0 + 7y v (k - Q) 


'K nT 


(£ f vf) 4 \fn 4(& f vf) 3 l e t]e 


log 


4 Ak\ 


(i s n ) + 0 (sb. (82) 


This gives Eq. (42) in the text. 

As in the case of 3D-AFQCP, the term -2ibnT / a 1 in Eq. (|421 ) cancels out the third term in Eqs. 


m an d dp . 

7. Derivation of Eq. (43) in the text 


Near the 2D-AFQCP, instead of Eq. (24) in the text, the correction to ASSO coupling a a {{rj) is 
given by 


, , f v (k,id) 2 J 

<*af(?7) = ---- + 


7v(k - Q) 
where cos p = (k • k')/(|k||k'|), and 


2 j r 2 "#_L r 
c* Jo 2 n Jl 


d y 


1 


y+^f(y 2 - 1) 




(83) 


a* = rj + K(l - cos p). (84) 

Let us introduce a A 1 !) expressing the second term on the r.h.s. of Eq. (1831) . Since the leading singu¬ 
larity for a An) in 77 arises from the integration with respect to ip near cos p ~ I of the term 1 / ^Ja*, in 
the denominator of integration with respect to y can be approximated by putting a* —> 0 : 

^ dp _1_ r°°Hv — 

JO 2n yjij + K{\ - cos p) J 1 y 2 

2 j r 2n dp 1 


a An) 


A/ r 2 
c* Jo 


I 


(85) 


C* VK Jo 2jt ^1 + (fj/K) - cos p 

The integration with respect to p in Eq. ll85l) is given by 2F (n/2, I / JI + (fj)2K)^ /n J2 + (rj/K), 
where F(tt/2, x) = K(x) is the complete elliptic function of first kind. The approximate form of K(x) 
near x ^ 1 is given by 


1 16 

Therefore, in the limit of 77 —» 0, we obtain 

. , , V27 , 32F 

a An) - ~ —p= log ——. 

nC* VK V 

Substituting definitions, C* = v 3 C/A, K = 2k^v 2 , and 77 = v^rj/A, into Eq. (1871) . we obtain 

7 2m* A 


a An) - 


2k“Vp TrvpC 


log 


6AAk 2 

n 


( 86 ) 


(87) 


( 88 ) 


Finally, by shifting 77 to 77 + Aq 2 to take into account the degree of distance from the hot point, we 
obtain the first term in the parenthesis on the r.h.s. of Eq. (43) in the text. 


13(20] 
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8. Derivation of Eq. (44) in the text 


Near the 2D-FQCP, instead of Eq. (37) in the text, the correction to ASSO coupling af(? 7 ) is given 


by 


win) = 


fy(k,ifl 

7v( k ) 

r y (k, id) 

7v(k) 


+ 27 


r>2n 

Jo 


dip \Jl - 


cos If 


2tt c* Vo* 


I 


dv 


V«*(l-cos If) 2 

7+ -— 7Z -(7 - 1) 


+ 2JS2d¥tn), 


(89) 


where C* = v 2 C/A. Substituting Eq. (l84l) into the integrands of the second term on r.h.s. of Eq. 
S 2 dF(??) is expressed as 

" 2rI dtp y/l ~ COS ip 


S 2 <i\ : (r]) - 


r 

Jo 


2?r C*yjfi + K{ 1 - COS 9?) 


r 


dv- 


7 + 


^l-cos ip -\Jfj+K( 1-cos <p) 2 


n2 ' 


(7 2 -l) 


(90) 


First, we calculate 5 2dF(0) which is considerably simplified: 

~' ln dip 1 


•S 2dF(0) = 


= ^~r f dy f 

c*Vkj i Jo 


2n- 


7+ -jr(l - cos ip)(y 2 - 1) 


n2 ’ 


(91) 


where we have interchanged the order of the integrations with respect to ip and y. The integration with 
respect to ip is easily performed to give 


S 2dF(0) = 


i r 
c*Vk Ji 


dy 


i y + ^(y 2 - 1) 


,3/2 




3/2' 


(92) 


With the use of a function I{x), defined by Eq. (45) in the text, S 2dF(0) is expressed as 


•WO) = jl 


Vk 

c* 


(93) 


Next, let us discuss the lowest-order correction in tj, which will turn out to be of the order of 
0{ \[rj). To this end, first, we change the integration variable in Eq. (1901) from ip to x = cos(q>/2): 

If 00 . 1 


S2dF(v) = 


V 2 r 1 

nC* J-i 


dx- 


yjrj + 2K{\-x 2 ) 


r 


dy 


y + Vl - x 2 -y/77 + 2K(l - x 2 ) (y 2 - 1 ) 


(94) 


Let us define the integration with respect to y in Eq. (l94l ) as M{x, fj): i.e., 

1 


M(x, f]) = J' 


dy 


y + Vl - x 2 y]f) + 2A^( 1 - x 2 ) (y 2 - 1) 


l2 ' 


(95) 


Then, 


S 2d¥{vi) ~ S 2dF(0) - 


t:C*Vk 


/' 


dx 


1 


Vl +{fil2K)-x 2 


M(x, fj) - 


1 


vr 


-.M(x, 0) 


(96) 


where we have used the fact that M{x, rj) is an even function in x, and changed the interval of the 
integration with respect to x from - I < x < I to 0 < x < I, multiplying the result by a factor 2. The 
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expression (l96l ) is rearranged as 

S 2dp(jl) ~ S 2dF(0) - 


Letter 


nC* VK 


[I'- 

•r 


l 


l 


dx- 


Vl + (rj/lK) - x 2 Vl -x 
' i [M(x, rj) - M(x, 0)] 


M{x, rj) 


10 VT-x" 

Let us introduce representing the first term in the brace of Eq. (f97l ). By making partial integra¬ 

tion with respect to x, Y^\rj) is expressed as 

1 


(97) 


Y a \v) = 


sin 


-l 


Vl + 0?/ 2K) 

l 


- sin -1 1 


M(l,fj) 


-f 


dx 

is reduced to 


sin 


Vl + W2K) 


- sin 1 x 


In the limit of rj —» 0, Eq. 

y(1)(77) " " V 5 M( 1,0) ~fo dx 


1 


>1 


+ 0{if) 


2 xfT^-2K 

By straightforward calculations, M( 1,0) and dM(x, 0)/dx are reduced to 


dM{x, rj) 
dx 


dMjx, 0) 
dx 


and 


f°° 1 

M(1,0)= dy — = 1, 

J t y 1 

dM{x, 0) 8 x[K 


dx 


C* 


r 


dy 


y + 2M(l - x 2)(y2 - 1 ) 


l3 


(98) 


(99) 


( 100 ) 


( 101 ) 


Namely, up to the leading order in tj, Y (l \r]) is given by 



F (1) (?7) - ~ 2 ^ + 0(ri)- 


( 102 ) 


It is easy to see that the second term in the brace of Eq. ( 1971 ) gives only the term of the order of 0(7]) 
at most. Therefore, S 2 dF(? 7 ) - S 2dF(0), Eq. (I97T ), is given by 

2 


S 2dl ; (//) - S 2dF(0) - - 


nx[2C*K 

1 V2A 


yjf] + 0(7)) 

x/rj + 0(rj), 


K 7rCvp 


(103) 


* = Cv^/A and 7] = v^r 

Finally, adding Eq. (1931 ) and Eq. (11031) and using the definition K = af(ij), Eq. (l89l ). is given 


where we have used the definitions, C* = Cv^/A and 7] = v^rj/A, in deriving the last equality. 


by 


, , f v (k,id) J 
a\(jl) = -+ 


7v(k) 


Ik^vl 


21 


VK) 2 V2A _ ^ 

■Fj-y57^ +ow 


(104) 


The second term on the r.h.s. of Eq. (11041) is nothing but the first and second terms in the bracket on 
the r.h.s. of Eq. (44) in the text because Vk/C* = V2 m*A/C. 


15 ® 
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9. Derivation of Eq. (46) in the text for 2D-AFQCP 

Here we derive Eq. (46) in the text for 2D-AFQCP. Instead of Eq. (fl 8 l) . r v (k, ie n ) in 2D is given 


by 


r v (k,ie n ) = Jy v (k - Q) 


nT 


Jo 2 /r 


1 


1 


1 


+ 


n’> 0 

+ 2 i b 


Vd+ (V^+ + e «') 2 VdE (VdT + e „') 2 

11 11 


(105) 


VH7 (V^+ + e «') 3 VdE (+ e„') 3 

where d+ = a + C*| + e„ - e„'| with a given by Eq. (15) in the text as a = ( v^ij/A ) + 2k 2 v£(l - cos p), 
and b - kpHjl _ cos t P)- 

Corresponding to Eq. ( fl9l ). the contribution from s n > = s n in the summation over s n > in the first 
and the third terms in the brace of Eq. (11051) is given by 


f v (k,i£n) = /yv(k-Q);rr 


Jo 2 /r 


1 


1 


+ 2 i b- 


1 


1 


(106) 


Vfl(V^ + £ n) 2 V^(V^ + £ n) 3 . 

Expanding the expression in the bracket up to linear order in e„ as in Eq. (l20l) . Eq. (11061) is reduced 


to 


f v (k, ie n ) = /y v (k - Q) nT 


n2n 

JO 




1 


2e„ 


+ 2 i 


o 27r \ [77 + K{ 1 - cos 9 c )] 3 / 2 [77 + AJI - cos p)) 2 

K (1 - cos p) 


_ 6 i [E. J— cos ^ )E,i } (107) 

2 [77 + K(l - cos p)) 2 V 2 [77 + K(l - cos <^)] 5 ^ 2 J 


The integration with respect to p of the first term in Eq. (11071) is performed as follows: 


r>2n 

Jo 


dp 


2n [77 + ATT - cos (c >)] 3 / 2 
d C 2n dp 


=(- 2) ^ r 

dri Jo 


1 


2 d f' 2n dp 


r 

Jo 


2tT y/fj + K( 1 - COS p) 

1 


(108) 


VZ di) Jo 2 n yi + (fi/K) - cos p 
The integration with respect to p in the last line of Eq. (11081) is the same as that appeared in Eq. 
(1851) and is given by 2 F [n/l, 1 / -J 1 + ( rj/lK )) r-J2 + (f)/K), where F(n/2, x ) = K(x) is the complete 
elliptic function of first kind and its approximate form near x - 1 is given by Eq. (l 86 l) . Therefore, in 
the limit 77 —> 0 , 


r»2n 

Jo 


dp 


1 1 , 32 

- =* -log ——-. 

Jo 27r yi + ifj/K) - cos p 7 rV 2 

Substituting this asymptotic form into the last line of Eq. (11081) . we obtain 


r*2 7T 

Jo 


1 ^ 1 V 21 

2n [77 + K{1 - cos (o )] 3 '' 2 y/~K 71 1 


dp 


(109) 


( 110 ) 
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The integration with respect to tp of the third term in Eq. (11071) is performed as follows: 


r-2n 

Jo 


d if 


(1 - cos If) 


o 2n \f) + K(l - cos ^)] 5 / 2 

_2\_d_ r 2n 

VdKj 0 


d if 


1 


2n [fi + K(l - cos </>)] 3 / 2 
Substituting the relation (II101) into Eq. (11111) . we obtain 

' 2n d f (1 - cos ip) 1 V 2 1 


f 


2n [77 + K{ 1 - cos <^)] 5 ^ 2 K 3 / 2 3 n fj 

The integration with respect to if of the last term in Eq. (11071) is performed as follows: 
" 2n dif 1 K + fj 


f 


2n [77 + K{\ - cos ip )] 2 
Then, to the leading order in 1 /?) » 1, we obtain 

’ 2n d ip 1 


[(K + f]) 2 - K 2 ) 3 ' 2 ' 


/ 


K 


(111) 


( 112 ) 


(113) 


(114) 


2/r [77 + K(l - cos if )] 2 {2Kfj ) 3 ' 2 ' 

The integration with respect to if of the third term in the bracket of Eq. (11071) is performed as 
follows: 


r 2n df b [k r 2n 
Jo W " Vi Jo 


d ip (1 - cos ip) 

2n [77 + K{\ - cos if)) 2 


d 

r 2n df 1 

-J* 

d 

1 

dK, 

Jo 2n f] + A 4 T - cos ip) 

V 2 

dK 

yj(K + rfj 2 - K 2 _ 


dK y]2K?i 


1 1 
4 K yfq 


(115) 


Substituting Eqs. (II101) . (II121) . (II141) . and (II151) into Eq. (11071) . we obtain 


71 Sji 


f v (k, ie n ) - Jy v ( k - Q) ( - 

Vk\ n V 2 t 12 


+ 1 


n 


2 i e„ 


2VK Vk f1 


(116) 


Substituting 77 = v^rj/A and K = 2k 2 v 3 into this expression, we obtain 
fy (k, ie n ) =* Jy v { k - Q) 


{Ak 2 ) 3/2 n T 
- e„ — --: ttt + 1- 


K nT 


Akl T 


Ak\ T _ v-f 

(kpvp) 3 q " 2 (k F v F ) 4 t ] 3/2 4 (k F v F ) 3 y/77 " (k F v F ) 4 77 


- 1 £„ 


(117) 


Since 77 oc T /(- log T) in the case of 2D-AFQCP. 2 ’ the first and the fourth terms in Eq. (II171) diverge, 
in the limit T —> 0, as (- log T), while the second term diverges as (- log T ’) 3/2 / 7" 1 2 _ \y c remark here 
that the first term in Eq. (II171) should be included in Eq. (43) in the text because it gives the same 
divergence in the T dependence as the first term in Eq. (43) in the text at least on the hot line (i.e., 
q m = 0) in the limit T —> 0. 

The summation over s n > (T s n ) in Eq. (11051) is performed in accordance with the 3D case, and. 
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instead of Eq. (l42l) . we obtain 

T v (k, ie n ) - f v (k, ie n ) + 7y v (k - Q) f ^ f 

Jo 2 n J n - 


Letter 


de" 


7TT 


Va + C*s"{ Va + C*s" + s") 2 


+ 7 tT + 2i b\ -- + — 


7lT 


a 2 a 2 


+ 0(s 2 n ), (118) 


where a = fj + K( 1 - cos <p) and b = yjK/2( I - cos </>)• Then, instead of Eq. (1431) . the linear terms 
of r v (k, ie n ) in e„ is obtained as follows: 


Ty (k, ie n ) T v (k, i£ n ) U„=0 ~ Ty (k, ifin) Ty (k, i£n) le„=0 

C 2n dip 

+Jy v ( k-Q) ^ 


nT --s n + 2^-{\s n ) 
a A a- 


■ (H9) 


Finally, substituting Eqs. (II141) . Ill 171) . and (II151) into Eq. (II191) . we obtain 

a 


Ty (k, i£ n ) - Ty(k, i£ n ) | e „ =0 + /y v (k - Q) 


Ak- n T 


Ak\ 1 


(&fVf ) 4 rj 4(^ f v f ) 3 y/rj 


(i s n ) + 0(4). ( 120 ) 


Here, we note that the term obtained from the first term in the square bracket of Eq. (II191) cancels out 
the second term in the bracket of Eq. (II171) . so that the term proportional to £„ disappears but only that 
proportional to i£„ remains. As in the case of 3D-AFQCP, this result verifies the fact that the first line 
in the brace of Eq. (fT 8 l) . i.e., its real part, is an even function in £„. Equation (11201) gives Eq. (46) in 
the text. 

As in the case of 3d-AFQCP, the term -2\bnT / a 2 in Eq. (II181) cancels out the third term in Eq. 
(fTTfil) or (1X171) . 

10. Derivation of Eq. (46) in the text for 2D-FQCP 

Here, we derive Eq. (46) in the text for the case of 2D-FQCP. We start with Eq. (11051) . the same 
equation as that in the case of 2D-AFQCP except that C*, in the parameter a+ = a + C*\ ±s n - £„/1, is 
not constant but C* = C* / yjl - cos Lp with C* being a constant. Then, we obtain the same expression 
for f v (k, i£ n ) as Eq. (11061) . However, since Eq. (11061) does not include C*, the expression Eq. (11071 ) 
for f y (k, i£ n ) can be applied to the present case of 2D-FQCP: 

' 2n dip j 1 2 £„ 


T v (k,i£ n ) = Jy v {k)nT 


J 

Jo 


2n \ [77 + K( 1 - cos 9 s )] 3 / 2 [77 + K{\ - cos </>)] 2 


+2i - 


(1 - cos <p) 


- 6 i 


2 [77 + K{\ - cos i ^>)] 2 V 2 [77 + K(\ - cos ^)] 5 ^ 2 
Then, the expression Eq. (II171) for f v (k, i£ n ) is also valid in the case of 2D-FQCP: 


( 121 ) 


T v (k,i£ n ) - /y v (k-Q) 


(Ak 2 ) 3/2 nT 


Ak* T v -vp7 

- £„ —-TTT +1- 


' Ak l nT 


Akl T 


~ 1£„ 


(k F v F ) 3 r\ 2(k F v F ) 4 rj 312 4(k F v F ) 3 y/rj (k F v F ) 4 77 
Since 77 oc (-T log T) in the case of 2D-FQCP, 3 ' the first and the fourth terms in Eq. (II171) vanish, in 


( 122 ) 


ism 
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the limit T —» 0, as l/(- log T), while the second term diverges as 1 /[ Vr(- log T) 2 ^ 2 ]. On the other 
hand, the first term in Eq. (11221) vanishes as l/(- log T) in the limit T —> 0, which should be included 
in Eq. (44) in the text because it gives a shaper cusp in T dependence in the limit T —* 0 than the 
second term of Eq. (44) in the text. 

The summation over s n ' (+ s n ) in Eq. (11051) is performed in the same way as the case of 2D- 
AFQCP, and we obtain the same relation as Eq. ( II191) : 


T v (k, ie n ) Ty (k, i£ n ) U„=0 — Ty (k, i^n) Ty (k, l£n) le„=0 



(123) 


where a = rj + K( I - cos ip) and b = ^K/ 2(1 - cos p). Then, T v (k, ie n ) is given by the same expression 


as Eq. (11201) : 


Ty (k, ie n ) - Ty (k, ie n ) | £n=0 + /y v (k - Q) 



This gives Eq. (46) in the text. 


As in the case of 2d-AFQCP, the term -2ibnT/a 2 in Eq. (II181) cancels out the third term in Eq. 
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